Nonequilibrium plasmons and transport properties of a double-junction quantum wire 
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We study theoretically the current-voltage characteristics, shot noise, and full counting statistics 
of a quantum wire double barrier structure. We model each wire segment by a spinless Luttinger 
liquid. Within the sequential tunneling approach, we describe the system's dynamics using a master 
equation. We show that at finite bias the non-equilibrium distribution of plasmons in the central 
wire segment leads to increased average current, enhanced shot noise, and full counting statistics 
corresponding to a super-Poissonian process. These effects are particularly pronounced in the strong 
interaction regime, while in the non-interacting case we recover results obtained earlier using detailed 
balance arguments. 
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I. INTRODUCTION 

The recent discovery of novel one-dimensional (ID) 
conductors that do not follow Fermi liquid theory has 
inspired extensive research activities both in theory and 
experiment. The generic behavior of electrons in ID con- 
ductors is well described by the Luttinger liquid (LL) the- 
ory which is ageneralization of the Tomonaga-Luttinger 
(TL) model |l|, |2J, |3|]. In one-dimensional conductors, 
unlike their higher dimensional counterparts, electron- 
electron (e — e) Coulomb interaction is poorly screened 
. Consequently, the fermionic quasiparticle excitations 
that are characteristic of Fermi liquids become unsta- 
ble in ID conductors: instead, collective density fluctua- 
tions constitute the stable elementary excitations in LLs. 
Luttinger liquids are further characterized by power-law 
correlations with interaction-dependent exponents, and 
by separation of the spin and charge degrees of freedom. 
Power-law behaviors of the differential conductance have 
been observed for edge states in the fractional quantum 
Hall regime p| and metallic single-walled carbon nan- 
otubes (SWNTs) 0, 0. The spin-charge separation 
was also observed in organic Bechgaard salts 0. 

One-dimensional single-electron tunneling transistors 
(SETs) that exhibit power-laws characteristic of Lut- 
tinger liquids have been fabricated using semiconducting 
quantum wires 01 or metallic SWNTs . In an exper- 
iment using semiconducting quantum wires, Auslaender 
et al. 10] showed that the widths of resonant levels of a 
ID island embedded in an interacting ID wire decrease 
as a power law over a range of temperatures, in a quan- 
titative _agreement with the theoretical prediction by Fu- 
rusaki 12]. In an SWNT experiment, on the other hand, 
Postma et al. studied a quantum dot (QD), formed 
by adjacent defects in a long metallic SWNT, in a SET 
geometry. The conductance as a function of temperature 
was seen to deviate from the conventional predictions 
[L2I Il3| . To explain the unpredicted temperature depen- 
dence of the conductance at low temperatures, Postma 



et al. [TJ, followed by Thorwart et al. [T^, proposed a 
new transport mechanism, correlated sequential tunnel- 
ing (CST), in which additional quantum correlations due 
to Coulomb interactions across the barriers were consid- 
ered beyond the conventional (uncorrelated) sequential 
tunneling (UST) approach. 

The power-law exponent of the temperature depen- 
dence of conductance in the UST and CST approaches in 
the strong interaction regime has been studied by several 
groups. While the numerical approach using a dynamical 
quantum Monte Carlo method supports CST approach 
[l5| , the "leading-log" methods followed by the functional 
renormalization group appr oaches does not support the 
CST mechanism ]ll [l3lll [13 . 

Since the pioneering works by Kane and Fisher [l3| . 
many properties of the double barrier (DB) structure 
with the Luttinger liquid leads have been investigated 
for quantum dots with single resonant level |2(j and with 
many resonant levels [T^ . Later, even the quantum dot 
was descried using the Luttinger model. In this regime, 
new phenomena arise due to the interplay between in- 
teractions within the ID wire and the Coulomb blockade 
induced by the confinement of electrons in a small region. 
Since the elementary excitations in the system are plas- 
mons (charge density waves), the excitation spectrum of 
the QD is bosonic. 

Very recently, various transport properties in such 
systems have been studied by a number of groups. 
Within the sequential tunneling approach, Braggio et 
al. found power-law-type differential conductance with 
sharp peaks related to the activation of plasmons |21j . 
charge-spin separation manifested in the conductance 
peak positions p^ . and shot noise indicating Luttinger 
liquid correlations |2^. The authors consistently as- 
sumed fast relaxation of the plasmonic excitations in the 
quantum dot, implying that excitations created by one 
tunneling event do not influence subsequent tunneling 
events. We hereafter refer to this approach as "equilib- 
rium plasmons" . 
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The present authors have focused more on the prop- 
erties and the consequences of the non-equillibrium dis- 
tribution of plasmons in the QD, in the following "non- 
equilibrium plasmons" p3, |25| . This is the case when 
the plasmon excitations in the quantum dot redistribute 
only via the single-electron tunneling events through tun- 
nel barriers. We found that while the steady-state plas- 
mon distribution in the QD is highly non-equilibrium, 
the average electric current is only weakly affected by 
the non-equilibrium properties of plasmons |24|. On the 
other hand, the non-equilibrium plasmons do affect more 
sensitive measurements such as shot noise (SN) and full 
counting statistics (FCS). Both of these characteristics 
show non-Poissonian behavior even at low bias voltages: 
the shot noise is greatly enhanced above the Poissonian 
limit (super-Poissonian) and the enhancement is more 
severe in the strong interaction limit |25|. 

As an extension of our previous work "Sinn, we inves- 
tigate the consequence of the non-equilibrium plasmons 
on average current, shot noise, and full counting statis- 
tics of a 1 D-SET that consists of three Luttinger liquid 
segments, based on a master equation approach in the 
conventional sequential tunneling regime (ST). In Sec. ^ 



we introduce our model for a ID quantum wire SET and 
present the analytically obtained tunneling rates within 
the golden rule approximation. In the same section, we 
also introduce the master equation which is used to ob- 
tain all the results of this work and discuss the possi- 
ble experimental realizations. Sec. I11II is devoted to the 
discussion of the distribution of the non-equilibrium oc- 
cupation probabilities of plasmonic many-body states. 
We then proceed to investigate the consequence of the 
non-equilibrium plasmons in the context of the transport 
properties. We hrst consider average current in Sec. IIVI 
and then discuss shot noise in Sec. E] Finally, we inves- 
tigate full counting statistics in Sec. IVII We conclude in 
Sec. En| 



II. FORMALISM 

The electric transport of a double barrier structure in 
the (incoherent) sequential tunneling regime can be de- 
scribed by the master equation [l2l 129. [23 
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where P(N, {n}, t) is the probability that at time t there 
are TV (excess) electrons and {n} = (711,712, • • • , n m , • • • ) 
plasmon excitations (i.e. collective charge excitations), 
that is, n m plasmons in the mode m on the quantum 
dot. The transitions occur via single-electron tunnel- 
ing through the left (L) or right (R) junctions (see Fig. 

The total transition rates T in master equation 
are sums of the two transition rates Tl and Tr where 
^l/r{N, {n} <— N',{n'}) is the transition rate from 
a quantum state (N', {n'}) to another quantum state 
(N, {n} via electron tunneling through L/R-] unction. 



Master equation Q implies that, with known tran- 
sition rates, the occupation probabilities P(N,{n},t) 
can be obtained by solving a set of linear first order 
differential equations with the probability conservation 
J2n {n} {n}, t) = 1. In the long time limit the sys- 
tem converges to a steady-state with probability distribu- 
tion limj-Kx, P(N, {n}, t + r) = P s t(N 7 {n}), irrespective 
of the initial preparation of the system. 



To calculate the transition rates, we start from the 
Hamiltonian of the system. The reservoir temperature is 
assumed zero (T = 0), unless it is stated explicitly. 



A. Model and Hamiltonian 

The system we consider is a (ID) quantum wire SET. 
Schematic description of the system is that a finite wire 
segment, which we call a quantum dot, is weakly coupled 
to two long wires as depicted in Fig. 2] The chemical 
potential of the quantum dot is controlled by the gate 
voltage (Vg) via a capacitively coupled gate electrode. In 
the low-energy regime, physical properties of the metal- 
lic conductors are well described by linearized dispersion 
relations near the Fermi points, which allows us to adopt 
the Tomonaga-Luttinger Hamiltonian for each wire seg- 
ment. We model the system with two semi-infinite LL 
leads and a finite LL for the central segment. The leads 
are adiabatically connected to reservoirs which keep them 
in internal equilibria. The chemical potentials of the 
leads are controlled by source-drain voltage (V) , and the 
wires are weakly coupled so that the single-electron tun- 
neling is the dominant charge transport mechanism, i.e. 
we are interested in the sequential tunneling regime. Rig- 
orously speaking, the voltage drop between the two leads 
(V) deviates from the voltage drop between the left and 
right reservoirs (say U) if electron transport is activated 
[28l l29j|. However, as long as the tunneling amplitudes 
through the junctions (barriers) are weak so that the 
Fermi golden rule approach is appropriate, we estimate 
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FIG. 1: (color online) Model system. Two long wires are adi- 
abatically connected to reservoirs and a short wire is weakly 
coupled to the two leads. Tunneling resistances at junction 
points Xl and Xr are Rl/r, an d the junction capacitances 
are considered equal Cl = Cr. Quantum dot is capacitively 
coupled to the gate electrode. 

The total Hamiltonian of the system is then given 
by the sum of the bosonized LL Hamiltonian Hq = 
Hl+Hd+Hr accounting for three isolated wire segments 
labeled by £ = (L, D, R), and the tunneling Hamiltonian 
Ht accounting for single-electron hops through the junc- 
tions L and R at and Xr respectively, 



H — Hq + H T 



(2) 



Using standard bosonization technique, the Hamilto- 
nian describing each wire segment can be expressed in 
terms of creation and annihilation operators for collec- 
tive excitations (St and b) HH3. For the semi-infinite 
leads, it reads 



M 



*P E 



for £ = L,R 



(3) 



where the index v labels the M transport sectors of the 
conductor and m the wave-like collective excitations on 
each transport sector. The effects of the Coulomb in- 
teraction in ID wire are characterized by the Luttinger 
parameter g v : g = 1 for noninteracting Fermi gas and 
< g < 1 for the repulsive interactions (g <C 1 in the 
strong interaction limit). The Coulomb interaction also 
renormalizes the Fermi velocity to v v = VFjg v . The en- 
ergy of an elementary excitation in sector v is given by 
e v = nhv u / L where L is the length of the wire and h the 
Planck constant. For instance, if the wire has a single 
transport channel (usually referred to as spinless elec- 
trons), e.g. a wire with one transport channel under a 
strong magnetic field, the system's dynamics is deter- 
mined by collective charge excitations (plasmons) alone 
(y = p and M = 1). If, however, the spin degrees of free- 
dom survive, the wire has two transport sectors (M = 2); 
plasmons (y = p) and spin- waves [y = a) [30| . If the sys- 
tem has two transport channels with electrons carrying 



spin (M = 4), as is the case with SWNTs, the trans- 
port sectors are total-charge-plasmons (u = p), rclative- 
charge-plasmons (y = Ap), total-spin- waves (y = a), and 
relative-spin- waves {y — Act) |3lLl33|. 

For the short central segment, the zero- mode need to 
be accounted for as well, which yields 



H 



D 



m. — l 



+ 2m { *'- NGl '' + ^m: N '-- E '- (4) 

In the second line of Eq. J2J, which represents the 
zero-mode energy of the quantum dot, the operator N v 
measures the ground state charge, i.e. with no excita- 
tions, in the ^-sector. The zero-mode energy systemat- 
ically incorporates Coulomb interaction in terms of the 
Luttinger parameter g p in the QD. To refer to the zero- 
mode energy later in this paper, we define the "charging 
energy" Ec as the minimum energy cost to add an excess 
electron to the QD in the off-Coulomb blockade regime, 
i.e., 

E c - [E d (2,0)-E d (1,0)] Ng=1/2 = e p /g. (5) 

Note that this is twice the conventional definition. The 
charging energy vanishes in the noninteracting limit 
(fjp = 1) an d becomes the governing energy scale in the 
strong interaction limit (g p <C 1). The origin of the 
charging energy in conventional quantum dots is the long 
range nature of the Coulomb interaction. In the theory 
of Luttinger liquid, the long range interaction can easily 
be incorporated microscopically through the interaction 
strength g p . For the effect of the finite-range interaction 
across a tunneling junction, for instance see Refs. I33T.I341 
Note that charge and spin are decoupled in Luttinger 
liquids, which implies the electric forces affect the (to- 
tal) charge sector only; due to intrinsic e — e interac- 
tions, g p < 1 but g U jtp = 1, and the gate voltage shifts 
the band bottom of the (total) charge sector as seen 
by the dimensionless gate voltage parameter Nq in Eq. 
0. As will be shown shortly, the transport properties 
of the L/R-leads are determined by the LL interaction 
parameter g p and the number of the transport sector 
M. In this work, we consider each wire segment has the 
same interaction strength for the (total) charge sector, 

g = g p L ^ — g p R ^ = gp D ^. Accordingly, the energy scales in 
the quantum dot are written by e p = — itTlvf / ' g P ^D 
and £o = e v^p = kTivf/Ld- 

We consider the ground state energy in the QD is the 
same as those in the leads, by choosing the reference 
energy E r in Eq. equals the minimum value of the 
zero-mode energy, 



. (e p N 2 G e p {l~N G ) 2 e (M-l) 
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\2Mg p ' 2Mo 



2M 



(6) 



where min(a;, x') denotes the smaller of x and x' , and the 
gate charge N G is in the range Nq € [0, M]. 
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The zero-mode energy in the QD, 



En — 



2Mg, 



■(N- 



^) 2 + ^E^ 2 -^ 



(7) 



yields degenerate ground states for N — and N = 1 
excess electrons when N G = [(M — l)g p + 1] /2. Here we 
replaced iV p by the number of the total excess electrons 
N since N p = N, and N v are all either even or odd 
integers, simultaneously; in the case of the SWNTs with 
N = ^2s.Ni t g excess electrons, where i — 1,2 is the 
channel index and s =|, | is the spin index of conduction 
electrons (M=4), N p = N,N a = '£ i {N i , t -N i , l ), N Ap = 



N, 



Es(Ni,s ~ N 2tS ), and N Ar7 = £ 4 (-l) 4 (A^ T 

From now on we consider only one spin-polarized (or 
spinless) channel unless otherwise stated — our focus is 
on the role of Coulomb interactions, and the additional 
channels only lead to more complicated excitation spec- 
tra without any qualitative change in the physics we ad- 
dress below. A physical realization of the single-channel 
case may be obtained e.g by exposing the quantum wire 
to a large magnetic field. 

For the system with high tunneling barriers, the elec- 
tron transport is determined by the bare electron hops at 
the tunneling barriers. The tunneling events in the DB 
structure are described by the Hamiltonian 



-L,R 



(8) 



where ^j(Xg) and &i(Xi) are the electron creation and 
annihilation operators at the edges of the wires near the 
junctions at Xl and Xr. As mentioned earlier, the elec- 
tron field operators & and W are related to the plasmon 
creation and annihilation operators b and b^ by the stan- 
dard bosonization technique. Different boundary condi- 
tions yield different relations between electron field oper- 
ators and plasmon operators. Exact solutions for the pe- 
riodic boundary condition have been known for decades 
0, l30j but the open boundary conditions which are apt 
for our system of consideration has been investigated only 
recently (See for example Refs. l3ll I35ll3ri 137^1 . 

The dc bias voltage V — Vl + Vr between L and R 
leads is incorporated into the phase factor of the tunnel- 
ing matrix elements ti = \ti\ exp^ieVit/h) by a time- 
dependent unitary transformation |38| . Here VlIr — 
VC/C R / L is voltage drop across the L/R-junction where 
C = ClCr/(Cl + Cr) is the effective total capacitance of 
the double junction, and the bare tunneling matrix am- 
plitudes \ti,m\ are assumed to be energy independent. 
Experimentally, the tunneling matrix amplitude is sen- 
sitive to the junction properties while the capacitance is 
not. For simplicity, the capacitances are thus assumed 
to be symmetric Cl = Cr throughout this work. By 
junction asymmetry we mean the asymmetry in (bare) 
squared tunneling amplitudes |£l/k| 2 . The parameter 
R = \~th\ 2 /\tFt\ 2 is used to describe junction asymmetry; 
R = 1 for symmetric junctions and R ^> 1 for a highly 
asymmetric junctions. 

It is known that, at low energy scales in the quantum 
wires with the electron density away from half-filling, the 



processes of backward and Umklapp scattering, whose 
processes generate momentum transfer across the Fermi 
sea (w 2kp), can be safely ignored in the middle of ideal 
ID conductors 0, including armchair SWNTs [3l|. The 
Hamiltonian J5J does not include the backward and Umk- 
lapp scattering (except at the tunneling barriers) and 
therefore it is valid away from half-electron-filling. 

We find that, in the regime where electron spin does 
not play a role, the addition of a transport channel does 
not change essential physics present in a single transport 
channel. Therefore, we primarily focus attention to a 
QW of single transport channel with spinless electrons 
and will comment on the effects due to multiple channel 
generalization, if needed. 



B. Electron transition rates 

The occupation probability of the quantum states in 
the SET system changes via electron tunneling events 
across L/R-j unctions. In the single-electron tunneling 
regime, the bare tunneling amplitudes I^l/aI are small 
compared to the characteristic energy scales of the system 
and the electron tunneling is the source of small pertur- 
bation of three isolated LLs. In this regime, we calculate 
transition rates T L / R between eigenstates of the unper- 
turbed Hamiltonian Hq to the lowest non- vanishing order 
in the tunneling amplitudes |£l/_r|. In this golden rule 
approximation, we integrate out lead degrees of freedom 
since the leads are in internal equilibria, and the transi- 
tion rates are given as a function of the state variables 
and the energies of the QD only [2ll l2~i| , 

r L /R (N',{n'}<-N,{n}) 

= Y\tL,R?l{W L/R )l D {{n>},{n}). (9) 

In Eq. © W L / R is the change in the Gibbs free energy 
associated with the tunneling across the L/R junction, 

W L/R = E D (N', {«'}) - E D (N, {n}) 

T (N> ~N)eV L/R . (10) 

Here E D (N,{n}) = (N,{n}\H D \N,{n}) is the energy of 
the eigenstate \N, {n}) of the dot and L/R correspond to 
— /+. For the QD with only one transport channel with 
spinless electrons only, 



E D (N, {n}) 



m—l 



(N - N G ) 2 
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D 



E r , (11) 



where e p — e p accounting that we consider only charge 



plasmons and n n 



M n b 



m I m m I "'TTl, 



is the number of plas- 



mons in the mode m. 

The function 7(2;) in Eq. 10 is responsible for the 
plasmon excitations on the leads, and given by (see e.g. 
Ref. ITT 



5 



7(£) = 2^ 



1 1 

2lT?i TTVp 



2^A g 1 /( 1 - 
Tlvf/3 



V 



+ a (3e 



-Pe/2 



r(l + a) ' 



(12) 



where /3 = 1/ksT is the inverse temperature in the leads, 
A is a short wavelength cutoff, and T(z) is the Gamma 
function. The exponent a = (g^ 1 — 1)/M is a character- 
istic power law exponent indicating interaction strength 
of the leads with M transport sectors (hence, in our case, 
M = 1). At g = 1 (non-interacting case), the exponent 
a = and it grows as g — ► (strong interaction) . The de- 
crease of the exponent a with increasing M implies that 
the effective interaction decreases due to multi-channel 
effect, and the Luttinger liquid eventually crosses over 
to a Fermi liquid in the many transport channel limit 

MM- 

For the non-interacting electron gas, the spectral den- 
sity 7(e) is the TDOS multiplied by the Fermi-Dirac dis- 
tribution function /fd(s) = [1 + exp(/3e)] _1 ; 7(e) = 
^K^fFDis) for g = 1. At zero temperature, 7(e) is 
proportional to a power of energy, 

Hm7( £ ) = e ( - £) ^-M^ 
t^o irhvF T(a 



1) 



(13) 



where ea — ^wf/A^ 1 ^ 1-9 - 1 is a high energy cut-off. At 
zero temperature 7(e) is the TDOS for the negative en- 
ergies and zero otherwise (as it should be), imposed by 
the unit step function 0(— e). 

The function j D in Eq. |JSJ accounts for the plasmon 
transition amplitudes in the QD, and is given by 

lD ({n'}, {n}) = S N ,, N+1 \(N', {n'}\*UXt)W, M)\ 2 
+ 5 N ,^ 1 \(N',{n'}\y D (X e )\N,{n})\ 2 (14) 



and L b a (x) are Laguerre polynomials. Additional trans- 
port sectors would appear as multiplicative factors of the 
same form as <f>({n?}, {n}) and result in a reduction of 
the exponent old- Notice that in the low energy scale 
only the first few occupations n m and n' m in the product 
differ from zero, participating to the transition rate © 
with nontrivial contributions. 



C. Plasmon relaxation process in the quantum dot 

In general, plasmons on the dot are excited by tun- 
neling events and have a highly non-equilibrium distri- 
bution. The coupling of the system to the environment 
such as external circuit or background charge in the sub- 
strate leads to relaxation towards the equilibrium. While 
the precise form of the relaxation rate, T p , depends on the 
details of the relaxation mechanism, the physical prop- 
erties of our concern do not depend on the details. Here 
we take a phenomenological model where the plasmons 
are coupled to a bath of harmonic oscillators by 



plasmon— bath 



EE 



m.n a 



,9mn b m b n a a 



+ h.c.) . (17) 



In Eq. (|17|) a a and aj, are bosonic operators describing 
the oscillator bath and g^ n is the coupling constants. We 
will assume an Ohmic form of the bath spectral density 
function 



J mn {uj) = Is™! 2 S{uj - u a ) = 7 p r w , 



(18) 



where we used that the zero-mode overlap is unity for 
N' = N±l and vanishes otherwise. The overlap integrals 
between plasmon modes are, although straightforward, 
quite tedious to calculate, and we refer to AppcndixIXlfor 
the details. The resulting overlap of the plasmon states 
can be written as a function of the mode occupations n m , 

K{n'}|$to)|{»})| 2 = \{{n'}\4>n(X e )\{n})\ 2 



L 



D 



^{n'},{n}) (15) 



with 
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(16) 



where oj a is the frequency of the oscillator correspond- 
ing to a a , 7 P is a dimensionless constant characterizing 
the bath spectral density, and T^" 1 = h 2 VFLD{\th\~ 2 + 
|iij|~ 2 ) is the natural time scale of the system. Within 
the rotating-wave approximation, the plasmon transition 
rate due to the harmonic oscillator bath is given by 



T p ({n'} <- {n}) = 7 p r 



Wp/e p 



(19) 



with W p = s P J2 m ( n 'm ~ n m), where e p = e p is the 
plasmon energy. Note that these phenomenological rates 
obey detailed balance and, therefore, at low temperatures 
only processes that reduce the total plasmon energy oc- 
cur with appreciable rates. 



D. Matrix formulation 



where n 



old 



«) 



(9i 



mhi(n' m ,n m ) and n 



(» 



max(nj n , n„ 



1) for the QD with one transport sector, 



For later convenience, we introduce a matrix notation 
for the transition rates T, with the matrix elements de- 
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fined by 



if (AO 



{«'}■{«} 



= Tt(N±l,{n'}*-N,{n}), (20) 



i.e, the element ({n 1 }, {n}) of the matrix block T e (N) is 
the transition rate Tg(N ± 1, {n'} *— N, {n}). Similarly, 



{n'},{n} 



Si n > 



{n'},{n} 



{„»} 



{n"},{n} 



. (21) 



and 



r P (A0 



= -r p ({n'},{n}) 

+ V}-WEW1'W)' ( 22 ) 

{»"} 



Master equation Q can now be conveniently expressed 
as 



-\P(t)) = -T\P(t)) 



(23) 



with r = r, 



■ p -r Ei=lm (P° - r t - f t ) , where \P(t)) 
is the column vector (not to be confused with the 
"ket" in quantum mechanics) with elements given by 
(N, {n} \P(t)) = P(N, {n} , t). Therefore, the time evo- 
lution of the probability vector satisfies 



\P(t)) = exp(-ri) |P(0)> . 



(24) 



In the long time limit, the system reaches a steady state 
|P(oo)). 

The ensemble averages of the matrices T e can then be 
defined by 



In contrast, the non-equilibrium occupation probabil- 
ity in the strong interaction regime with (nearly) sym- 
metric tunneling barriers depends only on the total en- 
ergy of the states, and follows a universal form irrespec- 
tive of electron charge N in the QD. In the leading order 
approximation, it is given by 



P<°»( £ ) 



1 



: exp 



3(a + 1) £ e 
log — 

2 eV B s n 



(26) 



where Z is a normalization constant. Notice that e is 
the total energy, including zero-mode and plasmon con- 
tributions. The distribution has a universal form which 
depends on the bias voltage and the interaction strength. 
The detailed derivation is in appendix [BJ This analytic 
form is valid for the not too low energies (e, eV) > 3e p 
and in the strongly interacting regime a > 1. More ac- 
curate approximation formula (Eq. i|B13(l ) is derived in 
appendix 151 

For symmetric junctions, the occupation probabilities 
fall on a single curve, well approximated by the analytic 
formulas Eqs. I|26(l and (|B 1.3(1 . as seen in the insets in 
Fig. [21 where P(s) is depicted as a function of the state 
energies for (a) g = 0.2 and (b) g = 0.5, with parameters 
R = 1 for the inset and R — 100 for the main figures 
(eV = 6e p , Nq = 1/2). For the asymmetric junctions, 
the line splits into several branches, one for each electric 
charge N, see the figure. However, as seen in[5[a), if the 
interaction is strong enough (g < 0.3 for R = 100 and 
eV = 6e p ), each branch is, independently, well described 
by Eq. I|2t)|) or (|B13|) . For weaker interactions, g > 0.3 
for R — 100 and eV = 6e p , the analytic approximation 
is considerably less accurate as shown in[2Ib). Even in 
the case of weaker interactions, however, the logarithms 
of the plasmon occupation probabilities continue to be 
nearly linear in e but with a slope that deviates from 
that seen for symmetric junctions. 



M*))= £ <iV,{n}|f, \P{t)) . (25) 

N,{n} 



We will construct other statistical quantities such as av- 
erage current and noise power density based on Eq. (|25|1 . 



III. STEADY-STATE PROBABILITY 
DISTRIBUTION OF NONEQUILIBRIUM 
PLASMONS 

By solving the master equation (|23l) numerically (with- 
out plasmon relaxation), in Ref. |2J, we obtained the 
occupation probabilities of the plasmonic many-body ex- 
citations as a function of the bias voltage and the in- 
teraction strength. We found that in the weak to non- 
interacting regime, a i=s or g = 1/(1 + a) ~ 1 for 
the wire with one transport sector, the non-equilibrium 
probability of plasmon excitations is a complicated func- 
tion of the detailed configuration of state occupations 
{n} = (ni,n 2 ,--- ,n m ,---). 



IV. AVERAGE CURRENT 

In terms of the tunneling current matrices Il/r across 
the junction L/R 



(27) 



the average current lL/Ft{t) = (Iz,/fl(*)) through L/R- 
j unction is 

h/ R {t)= (NAn}\lL/ R \P(t)) . (28) 

N,{n} 



Il/r - Te [Tl/r - t l/r) , 



The total external current I(t) = ^I(i)y, which includes 

the displacement currents associated with charging and 
discharging the capacitors at the left and right tunnel 
junctions, is then conveniently written as 



/(*) = E iPic t )i t {t), 



(29) 



e=L, R 
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FIG. 2: (color online) The occupation probability P(e) as a 
function of the mode energy e/e p . The energy e is abbre- 
viation for Ed(N, {n}), the bias eV = 6e p , the asymmetry 
parameter R = 100, and N G = 1/2 (T = 0). T wo a nalytic 
approximations, Eq. I2bll (blue curve) and Eq. 1B13II (cyan 
curve), are fitted to the probability distribution of the charge 
mode N — 1 (blue circle). The interaction parameter is (a) 
g = 0.2 and (b) g = 0.5. In the inset the case of symmetric 
junctions (R = 1) is plotted with the same conditions. 



where C^ 1 = 1 + C^ 1 . As the system reaches steady- 
state in the long time limit, the charge current is con- 
served throughout the system, / = /(oo) = Il(oo) — 

ifl(oo). 

One consequence of non-equilibrium plasmons is the 
increase in current as shown in Fig. [3| where the average 
current is shown as a function of the bias for different 
interaction strengths. The currents are normalized by 
I c = I(eV = 2Ec) with no plasmon relaxation (7^ = 0) 
for each interaction strength <?, and we see that the cur- 
rent enhancement is substantial in the strong interaction 
regime (g < 0.5), while there is effectively no enhance- 
ment in noninteracting limit g = 1 (the two black lines 
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FIG. 3: (color online) Average current I/I c as a function 
of the bias voltage eV/2Ec and LL interaction parameter g 
(black line for g = 1, green g = 0.7, red g = 0.5, blue g = 0.3, 
respectively) for R = 100 (highly asymmetric junctions) with 
no plasmon relaxation (7^ = 0, solid lines) or with fast plas- 
mon relaxation (70 = 10 4 , dashed lines). The bias voltage is 
normalized by the charging energy 2Ec and current is normal- 
ized by the current at eV = 2Ec with no plasmon relaxation 
for each g. Other parameters are Nq = 1/2, T = 0. 



are indistinguishable in the figure). In the weak interac- 
tion limit the current increases in discrete steps as new 
transport channels become energetically allowed, while 
at stronger interactions the steps are smeared to power 
laws with exponents that depend on the number of the 
plasmon states involved in the transport processes. 

Including the spin sector results in additional peaks 
in the average current voltage characteristic that can be 
controlled by the transverse magnetic field [22I lllj . 

The current-voltage characteristics show that, in the 
non-interaction limit, the non-equilibrium approach pre- 
dicts similar behavior for the average current as the de- 
tailed balance approach which assumes thermal equilib- 
rium in the QD. In contrast, in the strong interaction 
regime, non-equilibrium effects give rise to an enhance- 
ment of the particle current. 

Experimentally, however, the current enhancement 
may be difficult to attribute to plasmon distribution as 
the current levels depend on barrier transparencies and 
plasmon relaxation rates, and neither of them can be eas- 
ily tuned. We now turn to another experimental probe, 
the shot noise, which is more sensitive to non-equilibrium 
effects. 



V. CURRENT NOISE 

Noise in electronic conductors is given by the ensem- 
ble average of the current-current correlations |42l |43| . 
Thermal fluctuations and the discrete nature of the elec- 
tron charge are two fundamental sources of the noise; in 
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specific devices there may additional noise sources due 
to, e.g., fluctuating environmental variables. Thermal 
(equilibrium) noise is not very informative since it does 
not provide more information than the equilibrium con- 
ductance of the system. In contrast, shot noise, which 
is a consequence of the discreteness of charge and the 
stochastic nature of transport, can provide further insight 
beyond average current since it is a sensitive function of 
the correlation mechanism, internal excitations, and the 
statistics of the charge carriers 0, . 

Influence of quantum coherence on shot noise is an 
intriguing issue. It is known that the ensemble aver- 
aged quantum mechanical calculations of shot noise to 
the leading order is identical to the semiclassical ap- 
proaches when the comparable theories are available |45l 
in Sec. 5]. However, the word "semiclassical" should not 
be confused with the deterministic motion of the trans- 
port charges. For instance, Oberholzer et al. [4(| discuss 
the crossover from full quantum to classical shot noise, by 
tuning the electron dwell time in chaotic cavities, where 
by classical it means the deterministic nature of electron 
motion. 

Shot noise in interacting one-dimensional systems has 
also been the subject of many recent works. For instance, 
the shot noise of the edge states was used to measure the 
factional charge ge of the quasiparticles in the fractional 
quantum Hall states lil lia l50| . 

The shot noise of double-barrier structures was widely 
studied in the last decade. In conventional SET struc- 
tures, shot noise is suppressed below the Poisson limit 
due to the Coulomb correlations (in addition to the Fermi 
correlation): the e — e correlations typically result in 
reduction of shot noise. Both quantum mechanical ap- 
proaches |5lT l52| and semiclassical derivations based on 
a master e quat ion approach predict identical shot noise 
results |s3> I54L l55| , implying that the shot noise is not 
sensitive to the quantum coherence in DB structures. If 
the leads are superconducting in the SET structure, the 
tunneling particles are either single-electrons or Cooper- 
pairs and the shot noise is a functional of the dephasing 
process of the Cooper-pairs j5|| H3 • 

Semiclassical theories of shot noise based on a master 
equation approach in the sequential tunneling reg ime for 
a SET have been developed by many authors |53l |5J, |59j, 
I60l . l6ll.l62j . The predictions of some of these theories £3 
have been experimentally confirmed |63|. 

The noise power density in the steady state is given by 



be written as 59. 6 



S(oj) = lim 2 

t— »oo 



(It 



-\-10JT 



I(i + r)I(f))-(I(t) 



(30) 

The correlation functions Kw(t) = 
lim t _ s . 00 (ii(t + r)Ii'(t)^ can be deduced from the 
master equation (|23[) . In the matrix notation they can 



e 2 WWI 
N,{n} 

+ e(-T)U> ex P (+rv)i 



e(+r)i,exp(-rr)i^ 



|P(oo)> , (31) 



where 0(x) is the unit step function. 

To investigate the correlation effects, the noise power 
customarily compared to the Poisson value Sp i SSOn — 
2el. The Fano factor is defined as the ratio of the actual 
noise power and the Poisson value, 



2el ' 



(32) 



Since thermal noise (S = Ak B TG(V = 0)) is not 
particularly interesting, we focus on the zero frequency 
shot noise, in the low bias voltage regime where the 
Coulomb blockade governs the electric transport; T = 
and eV < 2E C . 

We begin by considering analytically tractable cases 
with only a few involved states, and then proceed to the 
full numerical results. The finite frequency shot noise is 
briefly discussed in subsection IV El 



A. Two-state model; eV « e p 

The electron transport involving only two lowest en- 
ergy states in the quantum dot are well studied by many 
authors (see, for instance, Ref. liBj) . Nevertheless, for 
later reference we begin the discussion of shot noise with 
two-state process, which provides a reasonable approx- 



imation for eV 



e„. At biases such that eV 



and sufficiently low temperatures, the two lowest states 
I A, ni) = |0, 0) and |1, 0) dominate the transport process 
and the rate matrix is given by 



- 7 + 7 - 



(33) 



where the matrix elements are 7+ = Tl(1, *— 0, 0) and 

7- = r fl (o,o«-i,o). 

With the current matrices defined by Eq. I|27|l 





7+ 



7- 





the noise power is obtained straightforwardly by Eqs. 
(|3U|I and (|31l) using the steady state probability 



|P(oo)) = 



■Poo 
Pud 



1 



7" 



The Fano factor (|32|) takes a simple form 

1 + (7~/7 + ) 2 
(1 + 7-/7+) 2 



F 2 = R 



no 



P?o 



(34) 



(35) 
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where 



is 



7- _ 1 (eV/2-5E 
7+ ~ R \eV/2 + SE 

and SEq is the shift of the bottom of the zero mode energy 
induced by the gate voltage, 



SE Q = {8N G )e p /g, 5N G = N G - 1/2. 



(36) 



Note that Eq. is valid for \5E \ < eV/2, other- 

wise 1 = and S = due to Coulomb blockade. We 
see from Eq. (|35fl that the Fano factor is minimized for 
7~/7 + = 1 and maximized for (7~/7 + ) ±:L = 0, with the 
bounds 1/2 < F 2 < 1. At the gate charge N G = 1/2, 
it is determined only by the junction asymmetry param- 
eter R: F 2 = (1 + i? 2 j/(l + i?) 2 . Note that when only 
two states are involved in the current carrying process 
(ground state to ground state transitions), the Fano fac- 
tor cannot exceed the Poisson value F = 1. 

As a consequence of the power law dependence of the 
transition rates on the transfer energy O, the Fano 
factor is a function of the bias voltage, gate voltage, and 
the interaction strength; it varies between the minimum 
and maximum values 



1 

2' 
1, 



at 8E 



1 - R 1 ' 01 eV 



at 5E =±^- 



(37) 



(cf. Eq. (5) in Ref. |23|). Fig. gjdepicts the Fano factor as 
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FIG. 4: (color online) Fano factor F = S(0)/2eI as a function 
of the gate charge Ng and LL interaction parameter g (black 
line for g = 1, green g = 0.7, red g = 0.5, blue g = 0.3, ma- 
genta g = 0.2, respectively) for R = 100 (highly asymmetric 
junctions) at eV = e v (T = 0). 



a function of gate charge N G and interaction parameter 
g at the bias eV = e p for strongly asymmetric junctions 
R = 100 (T = 0). The gate charge N G corresponding to 
the minimum of the Fano factor (F = 1/2) in the figure 



N* lp = l/2 + g(eV/2e p ) 



1 - i? 1 /" 

1 + W7°> ' 



(38) 



The Fano factor independently of the interaction strength 
crosses F 2 = (1 + i? 2 )/(l + R) 2 at N G = 1/2, and it 
approaches maximum F 2 = 1 at 2Vg = 1/2 ± <7(ey/2e p ). 



B. Three-state model; eV ^ 2e p 

The two-state model is applicable for bias voltages be- 
low eVth = 2{e p —\8Eo\), since at least three states can be 
involved in transport above this threshold voltage. For 
electron transport involving three lowest energy states 
in the quantum dot, \N,ni) = |0,0),|1,0), and |1,1) 
with n m = for m > 2, the noise power can be calcu- 
lated exactly if the the contribution from the (backward) 
transitions against the bias is negligible, as is typically 
the case at zero temperature. In practice, however, the 
backward transitions are not completely blocked for the 
bias above the threshold voltage of the plasmon excita- 
tions, even at zero temperature: once the bias voltage 
reaches the threshold to initiate plasmon excitations, the 
high energy plasmons in the QD above the Fermi ener- 
gies of the leads are also partially populated, opening the 
possibility of backward transitions. 

A qualitatively new feature that can be studied in the 
three-state model as compared to the two-state model 
is plasmon relaxation: the system with a constant total 
charge may undergo transitions between different plas- 
mon configurations. 

We will show in this subsection that the analytic solu- 
tion of the Fano factor of the three-state process yields an 
excellent agreement with the low bias numerical results 
in the strong interaction regime, while it shows small 
discrepancy in the weak interaction regime (due to non- 
negligible contribution from the high energy plasmons). 
We will also show that within the three state model the 
Fano factor may exceed the Poisson value. 



1. Analytic results 

By allowing plasmon relaxation, the rate matrix in- 
volving three lowest energy states \N, n\) = |0, 0) , |1, 0) , 
and |1, 1) is given by 



7o 



+ 7i 

7o + 

it 



-To ~7i 
7o _-7p 
7i +7p 



(39) 



with the matrix elements 7^ = Fl(1,« <— 0,0), 7" = 
r fl (0, <— 1, i), i = 0, 1 and 7 P introduced in Eq. (T§Jl. 
Current matrices defined by Eq. (|27|l are 

[tt]{2,l} = To" . [it]{3,l} = Tl": [iil]{l,2} = To ! [itf]{l,3} = Tl~ 
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with = for other set of i,j — 1,2,3, where 

t = L,R. The noise power is obtained straightforwardly 
by Eqs. (|3T]|) and l|3TJ) using the steady state probability 



|PM) = 



-Poo 
Pio 

Pll 



1 

Z 



7o(7i +7p) 
7(1(71 +lp)_+ltl P 
7i + 7o 



(40) 



with normalization constant Z = 7^7^ +77/7(7 +7c/77/ + 
(70 + 7(7 + lt)lp- Using the average current / = 
e(-Pio7o~ + Pii7i~)> tn e Fano factor F = 5(0)/2eJ is given 
by 



^ =P, 



00 



P 
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p2 i + 2Pl 12b 

Z 7o 



+ 



(7o ) 2 + (7i ) 2 - 7o 7i + 



7o + 7i - 

+ — 7i 7p 

7o 



(41) 



Compared to the Fano factor l|35[l in the two-state pro- 
cess, complication arises already in the three-state pro- 
cess due to the last term in Eq. I|41[) which results from 
the coupling of P\\ and the rates which cannot be ex- 
pressed by the components of the probability vector. 

In order to have \N, rix) — |0, 0) , |1, 0) , and |1, 1) as the 
relevant states, we assume 70 > 7^ or more explicitly 
Ng > Ng p which is introduced in Eq. fifty • In the 
opposite situation (7^7 < 7(7), the relevant states are 
\N, m) = |0, 0) , |0, 1) , and |1, 0), and above description is 
still valid with the exchange of electron number N = <-> 
1 and the corresponding notations 7+ <-> 7^, i = 0, 1. 
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0.8 



FIG. 5: (color online) Fano factor F = S(0)/2eI as a func- 
tion of the gate charge Ng for symmetric junctions (R = 1) 
at voltage eV = 2e p (T — 0), with no plasmon relaxation. 
Numerical results (solid lines) vs. analytic results with three 
states, Eq. 1411 (dashed lines) for g — 0.3 (blue), 0.5 (red), 
and 1.0 (green). 

To see the implications of Eq. (|41|l , we plot the Fano 
factor in Fig. [S] with respect to the gate charge Ng 
for symmetric junctions at eV = 2e p (T = 0), with no 
plasmon relaxation (7 P = 0). 



Two main features are seen in Fig. [S] Firstly, the 
shot noise is enhanced over the Poisson limit (F = 1) 
in the strong interaction regime, g < 0.5, for a range of 
parameters with gate charges near (but not including) 
Ng = 1/2. As discussed above, in the low bias regime 
eV < 2(s p — \SEq\) at zero temperature, no plasmons are 
excited and the electric charges are transported via only 
the two-state process following the Fano factor (|35|l which 
results in the sub-Poissonian shot noise (1/2 < F < 1). 
Once the bias reaches the threshold eVth = 2(e p — |(W?o|), 
it initiates plasmon excitations which enhance the shot 
noise over the Poisson limit. This feature is discussed in 
more detail below. 

Secondly, in the weak interaction regime (g 0.5) 
a small discrepancy between the analytic result 141|) 
(dashed line) and the numerical result (solid line) is 
found. It results from the partially populated states 
of the high energy plasmons over the bias due to non- 
vanishing transition rates. On the other hand, a sim- 
ple three-state approximation shows excellent agreement 
in the strong interaction regime (<? = 0.3 in the fig- 
ure), indicating negligible contribution of the high energy 
plasmons (Ed(N, {n}) > eV/2) to the charge transport 
mechanism. This is due to the power law suppression 
of the transition rates l|13|) as a function of the transfer- 
energy (|Trj|) . 



2. Limiting cases 

To verify the role of non-equilibrium plasmons as the 
cause of the shot noise enhancement, we consider two 
limiting cases of Eq. (|4*T|l : 7 p = and j p ^> 7* . 

In the limit of no plasmon relaxation (7 P = 0), the 
Fano factor i|41|) of the three-state process is simplified 
as 



P 3 (0) -l-2[(l 



Po)P 



10 



2Pi Pi 



(l-P u )P n ] 

(7o~) 2 + (7r) 2 



7o 7i 



with steady-state probability 



|P(oo)> = 



Poo 
Pro 

Pll 



1 

Z 



7o 7i_ 
7(7 7i_ 
7^70 



(42) 



(43) 



where the new normalization constant is Z 

v + 



7o 7i 



7i 7 + 7 7i • 

The three-state approximation is most accurate in the 
low bias regime 2 (e p — \ SE 1) < eV <C 2E C (SE is defined 
in Eq. Ip^Hjl l and for gate voltages away from Nq = 1/2, 
i.e., for 1/2 < N G < 1 (or < N G < 1/2 with the 
exchange of indices regarding particle number N — 
1). In this regime, 7^ and/or 7.7 dominate over the other 
rates, (7^, 777) > (7(7,7i")> which results in 1 > Pi > 
(Po, P11) w 0. The Fano factor (|42J) now reduces to 



F. 



(0) 



l + 2P 10 Pn- 



(7q-) 2 + (7D S 
7(7 77/ 



1 



>1L 

'it 



(44) 
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Notice that while Eq. I|42|) is an exact solution for the 
three-state process with no plasmon relaxation, Eq. I)44|) 
is a good approximation only sufficiently far from Nq — 
1/2. In this range, Eq. (|44ll explicitly shows that the 
opening of new charge transport channels accompanied 
by the plasmon excitations causes the enhancement of the 
shot noise (over the Poisson limit). 

In the limit of fast plasmon relaxation, on the other 
hand, "f p 3> 7_- and effectively no plasmon is excited, 



|P(oo)> 



Poo 
Pio 

Pll 



7o + 



it 



7o 



7o 

7o + +7i + 




Consequently, the Fano factor (|41|l is given by 

2 



F, 



(oo) 



2 P ni 



e[il]. 



(45) 



(46) 



The maximum Fano factor F = 1 is reached if one of 
the rates 7^ or 7^ dominates, while the minimum value 
F = 1/2 requires that 7^ = 7o + 7i , i.e., that the 
total tunneling-in and tunneling-out rates are equal;more 
explicitly, 



1 

R 



1 - 



26E 



eV/2 + SE C 



1 - 



eV/2 + SE 



= 1. 



(47) 



Numerical results 



The limiting cases of no plasmon relaxation (7 P = 0) 
and a fast plasmon relaxation (7^ = 10 4 ) are summa- 
rized in Fig. Ela) and (b), respectively. In the figure the 
Fano factor is plotted as a function of gate voltage and 
interaction parameter g in the strong interaction regime 
0.6 < g < 0.3 at eV = 2e p for R = 100 (strongly asym- 
metric junctions). 

As shown in Fig. HJa), the Fano factor is enhanced 
above the Poisson limit (F = 1) for a range of gate 
charges away from Ng = 1/2, especially in the strong in- 
teraction regime, as expected from the three-state model. 
The shot noise enhancement is lost in the presence of a 
fast plasmon relaxation process, in agreement with ana- 
lytic arguments, as seen in Fig. Elb) when F is bounded 
by 1/2 < F < 1. Hence, slowly relaxing plasmon excita- 
tions enhance shot noise, and this enhancement is most 
pronounced in the strongly interacting regime. 

In the limit of fast plasmon relaxation F exhibits a 
minimum value F = 1/2 at positions consistent with 
predictions of the three-state model: the voltage polarity 
and ratio of tunneling matrix elements at the two junc- 
tions is such that total tunneling-in and tunneling-out 
rates are roughly equal for small values of Ng- If plas- 
mon relaxation is slow, F still has minima at approxi- 
mately same values of Ng but the minimal value of the 
Fano factor is considerably larger due to the presence of 
several transport channels. 




(a) 




(b) 

FIG. 6: (color online) Fano factor F = S(0)/2eI as a function 
of the gate charge Ng and LL interaction parameter g for 
R — 100 (highly asymmetric junctions) at eV = 2e p (T = 0), 
(a) with no plasmon relaxation (7^ = 0) and (b) with fast 
plasmon relaxation (7^ = 10 4 ). 



D. Interplay between charge fluctuations and 
plasmon excitations near eV = 2Ec 



So far, we have investigated the role of non-equilibrium 
plasmons as the cause of the shot noise enhancement and 
focused on a voltage range when only two charge states 
are significantly involved in transport. The question nat- 
urally follows what is the consequence of the charge fluc- 
tuations. Do they enhance shot noise, too? 

To answer this question, we first consider a toy model 
in which the plasmon excitations are absent during the 
single-charge transport. In the three- N-state regime 
where the relevant states are \N) = | — 1),|0) and |1) 
with no plasmon excitations at all. At zero temperature, 
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the rate matrix in this regime is given by 



rt! -r " o 
-rii r - + r+ -17 
-r+ rr 



(48) 



where the matrix elements are = + 1,{0} 
i,{0}),Tr = r R (i-l,{0} «- t,{0}), i = -1,0,1. Re- 
peating the procedure in subsection IV Bl we arrive at a 
Fano factor that has a similar form as Eq. (|42|) . 



= 1 - 2[(1 - P_i)P_i + (1 - P x )Pi] 



2P-iP 



, 1 1 , 

with the steady-state probability vector 



(49) 



|P(oo)) = 



P-x 
Po 

Pi 



1 

Z 



r -rr 

r± x rr 

r+ r+ 
1 -l 1 



(50) 



where Z = T^T^ + TtjT^ + T+ 1 T+. 

Despite the formal similarity of Eq. l|4*9"]) with Eq. l(l2*jl. 

its implication is quite different. In terms of the transi- 
tion rates, F^n reads 



F- 



3N 



2 



rrro(ri 1 r+ + (T± 1 -r+)rr) 
•riir+(rorr + (rr-ro)rr 



(51) 



Since Tt-^ > and T 1 > T Q in the three-N-state 
regime, the Fano factor F 3 jy is sub-Poissonian, i.e., 
F3N < 1, consistent with the conventional equilibrium 
descriptions [23ll59]|. 

We conclude that while plasmon excitations may en- 
hance the shot noise over the Poisson limit, charge fluctu- 
ations, in contrast, do not alter the sub-Poissonian nature 
of the Fano factor in the low energy regime. This quali- 
tative difference is due to the fact that certain transition 
rates between different charge states vanish identically 
(in the absence of co-tunneling): it is impossible for the 
system to move directly from a state with N = — 1 to 
TV = 1 or vice versa. 

Therefore, we expect that for bias voltages near eV — 
2Ec > 2e p , when both plasmon excitations and charge 
fluctuations are relevant, the Fano factor will exhibit 
complicated non-monotonic behavior. Exact solution is 
not tractable in this regime since it involves too many 
states. Instead, we calculate the shot noise numerically, 
with results depicted in Fig. [7| where the zero tempera- 
ture Fano factor is shown as a function of the bias eV and 
LL interaction parameter g for R = 100 at Nq — 1/2, (a) 
with no plasmon relaxation (7 P = 0) and (b) with fast 
plasmon relaxation (7^ = 10 4 ). 

In the bias regime up to the charging energy eV < 
2Eq, the Fano factor increases monotonically due to non- 
equilibrium plasmons. On the other hand, the charge 




eV/E c 2 



(a) 




(b) 

FIG. 7: (color online) Fano factor F = S(0)/2e7 as a function 
of the bias eV and LL interaction parameter g for R = 100 
(highly asymmetric junctions) at Ng = 1/2 (T — 0), (a) with 
no plasmon relaxation ("/ p = 0) and (b) with fast plasmon 
relaxation (7^ = 10 4 ). 



fluctuations contribute at eV > 2Eq which tends to sup- 
press the Fano factor. As a consequence of this competi- 
tion, the Fano factor reaches its peak at eV = 2Ec and 
is followed by a steep decrease as shown in Fig. E^a). 
Note the significant enhancement of the Fano factor in 
the strong interaction regime, which is due in part to the 
power law dependence of the transition rates with expo- 
nent a = (1/g — 1) as discussed earlier, and in part to 
more plasmon states being involved for smaller g since 
Ec = £p/g- The latter reason also accounts for the fact 
that the Fano factor begins to rise at a lower apparent 
bias for smaller g: the bias voltage is normalized by Ec 
so that plasmon excitation is possible for lower values of 
eV/Ec for stronger interactions. 

In the case of fast plasmon relaxation, the rich struc- 
ture of the Fano factor due to non-equilibrium plasmons 
is absent as shown in Fig. Efb), in agreement with the 
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discussions in previous subsections. The only remaining 
structure is a sharp dip around eV — 2E C that can be 
attributed to the charge fluctuations at eV > IEq. Not 
only is the minimum value of the Fano factor a function 
of the interaction strength but also the bias voltage at 
which it occurs depends on g. The minimum Fano fac- 
tor occurs at higher bias voltage, and the dip tends to 
be deeper with increasing interaction strength. Note for 
g < 0.5, the Fano factor did not reach its minimum still 
at largest voltages plotted (eV/(2E c ) = 1-2). 

The Fano factor at very low voltages eV < 2e p for 
N G = 1/2 is = (1 + R 2 )/{1 + R) 2 (sec Eq. ©), re- 
gardless of the plasmon relaxation mechanism. As shown 
in Fig. Etb), the Fano factor is bounded above by this 
value in the case of fast plasmon relaxation. Notice that 
in the case of no plasmon relaxation Qa)), the dips in 
the Fano factor at eV > 2Ec reach below . See Fig. 
1(a) in Ref. |25| for more detail. 

Since both the minimum value of F and the voltage at 
which it occurs are determined by a competition between 
charge fluctuations and plasmonic excitations, they can- 
not be accurately predicted by any of the simple analytic 
models discussed above. 



E. Finite frequency noise 

In the high frequency limit, lu — * oo, the correlation 
effects are lost in the noise power l|30|l except the S(t) - 
term in Eq. (|31l) that reflects the Pauli exclusion [6(J, 
and the asymptotic value of the noise spectrum reduces 
to 

(52) 

where 

A t = e £ (JV,{n}|fJ + f;|P(oo)) (53) 

N,{n} 

is the total tunneling rate across the junction I without 
regard to direction. 

The decay of the current-current correlations at finite 
frequencies is depicted in Fig. where the shot noise 
power S(uj)/2el is shown as a function of the bias and 
the frequency for a strong interaction (g — 0.5) and the 
strongly asymmetric tunnel barriers (R = 100) (a) with 
no plasmon relaxation and (b) with fast plasmon relax- 
ation (7,3 = 10 4 ). 

In the regime of the elastic process, in which the charge 
transport does not involve excitations (plasmons) in the 
dot, the backward tunneling against the bias is blocked 
at zero temperature. In this regime, the rate density is 
equal to the average current, Al = Ar — /, and the 
high frequency noise power asymptotically converges to 
its minimum value S(oo) = el. This is clearly seen in 
the regime of the low bias eV < 2e p (= Ec) in Fig. 02 a ) 
and (b). 




(b) 

FIG. 8: (color online) Finite frequency shot noise S(u>)/2el as 
a function of the bias eV/Ec and frequency ecu /I for R = 100 
and g = 0.5, at Nq = 1/2 (T = 0), (a) with no plasmon 
relaxation (-y p = 0) and (b) with fast plasmon relaxation 
(7p = 10 4 ). 



In the presence of fast plasmon relaxation, still no 
backward tunneling is possible and the shot noise reaches 
its minimum value S(oo) = el at high frequency limit 
(see Fig. Etb)). However, the non-equilibrium plasmons 
excited above the Fermi energies of the leads invoke non- 
vanishing backward tunneling and the high frequency 
shot noise remains above e/, as shown in the bias regime 
of eV > 2e p in Fig. 02a). 

Another identifiable feature in Fig. IBfa) is the rapid 
decrease of the shot noise power as a function of uj at 
those voltages when the Fano factor has a maximum. 
The maxima occur at voltages when additional states be- 
come significantly populated, and the characteristic low 
frequency reflects the slow transition rates to states near 
their energy thresholds. 
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VI. FULL COUNTING STATISTICS 



tion 



Since shot noise that is a current-current correlation 
is more informative than the average current, we expect 
more information with higher-order currents or charge 
correlations. The method of counting statistics, which 
was introduced to mesoscopic physics by Levitov and 
Lesovik [64| followed by Muzykantskii and Khmelnitskii 
[65j and Lee et al. [6(| , shows that all orders of charge cor- 
relation functions can be obtained as a function related to 
the probability distribution of transported electrons for 
a given time interval. This powerful approach is known 
as full counting statistics (FCS). The first experimental 
study of the third cumulant of the voltage fluctuations 
in a tunnel junction was carried out by Reulet et al. [6?| . 
The experiment indicates that the high cumulants are 
more sensitive to the coupling of the system to the elec- 
tromagnetic environment. See also Levitov and Lesovik 
[68| and Beenakker et al. [6!| for the theoretical discus- 
sions on the third cumulant in a tunnel-barrier. 

de Jong [7Cj showed that the low-energy physics of 
the FCS calculations both with quantum mechanical and 
semiclassical approaches are identical in the DB struc- 
tures, as in the explicit shot noise calculations, except 
for a short initial time scale. 

We will now carry out an FCS analysis of transport 
through a double barrier quantum wire system. The 
analysis will provide a more complete characterization 
of the transport properties of the system than either av- 
erage current or shot noise, and shed further light on the 
role of non-equilibrium vs. equilibrium plasmon distri- 
bution in this structure. 

Let P(M, t) be the probability that M electrons have 
tunneled across the right junction to the right lead during 
the time r. We note that 



P(M,t)= lim V 

t— »00 * » 

N,{n} M ,N ,{n } 



E 



P(M + M, N, {n},t + r; M , N , {n Q } , t), 

where P(M + M, N, {n} , t + r; M , N , {n } , t) is called 
joint probability since it is the probability that, up to time 
t, Mo electrons have passed across the right junction and 
Nq electrons are confined in the QD with {no} plasmon 
excitations, and that Mq + M electrons have passed It- 
junction with (JV, {n}) excitations in the QD up to time 
t + r. The master equation for the joint probability can 
easily be constructed from Eq. (|23|l by noting that M — > 
Mil as iV ^ TV =F 1 via only R-junction hopping. 

To obtain P(M,t), it is convenient to define the char- 
acteristic function conjugate to the joint probability as 

g(9 7 N 7 {n} 7 r) = lim £ e~" M ]T 

M M ,N ,{n } 

P(M + M, N, {n},t + r; M , N , {n Q } , t). 
The characteristic function satisfies the master equa- 



— \g(9,r)) = -T(9) \g(9,r)) 



(54) 



with the initial condition \g(9, r = 0)) = |P(oo)). The 
^-dependent T in Eq. I|54|l is related to the previously 
defined transition rate matrices through 



f(0) = f>(f° - rj - r~ 



(55) 

The characteristic function G(9, r) conjugate to P(M, r) 
is now given by G(9, r) = £ N ,{n} W \ n ] \9( 9 > T )> - or 



N,{n} 



T(9)t |P(oo)) . (56) 



Finally, the probability P(M, r) is obtained by 



P(M,r) 



dS 
2^ 



\-iOM 



G(9,t) 



dz G(z,t) 
2ri z 2A/ +! 



(57) 

with z = e l6 / 2 y where the contour runs counterclockwise 
along the unit circle and we have used the symmetry 
property G(z,t) = G(—z,r) for the second equality. 

Taylor expansion of the logarithm of the characteristic 
function in i9 defines the cumulants or irreducible corre- 
lators Kfe(r): 



lnG(fl,r)=f;^/ Sk (r) 



fe=i 



(58) 



The cumulants have a direct polynomial relation with 
the moments n k (r) = ~^2 n n k P(n, t). The first two cu- 
mulants are the mean and the variance, and the third 
cumulant characterizes the asymmetry (or skewness) of 
the P{M, t) distribution and is given by 



^3 (t) = Sn(r) 3 = (n(r) — n(r)) 2 



(59) 



In this section, we investigate FCS mainly in the con- 
text of the probability P(M, t) that M electrons have 
passed through the right junction during the time r. 
Since the average current and the shot noise are propor- 
tional to the average number of the tunneling electrons 
(M) and the width of the distribution of P(M, r), respec- 
tively, we focus on the new aspects that are not covered 
by the study of the average current or shot noise. 

In order to get FCS in general cases we integrate the 
master equation l|54|l numerically (see subsection IVI C|l . 
In the low-bias regime, however, some analytic argument 
can be made. We will show through the following subsec- 
tions that for symmetric junctions in the low-bias regime 
(2e p < eV < 2Ec), and irrespective of the junction sym- 
metry in the very low bias regime (eV < 2e p ), P(M,t) 
is given by the residue at z = alone, 



P(M,r) 



1 



i2M 



(2M)! dz 2M 



G(z,t). 



(60) 



15 



Through this section we assume that the gate charge is 
Nq = 1/2, unless stated explicitly not so. 

We will now follow the outline of the previous section 
and start by considering two analytically tractable cases 
before proceeding with the full numerical results. 



A. Two-state process; eV w e p 

For the very low bias eV < 2e p at zero temperature, no 
plasmons are excited and electrons are carried by tran- 
sitions between two states (N,n) = (0,0) <-> (1,0). In 
this simplest case, the rate matrix T(8) in Eq. 154(1 
is determined by only two participating transition rates 
7+ = T L (1, <- 0, 0) and 7 - = r«(o, <- 1,0), 



r(0) 



r Q + 5 -(T -8)e- 
-(T + S) T -5 



(61) 



with T = ( 7 + + 7 ")/2, S = ( 7 + - 7 ")/2. 

Substituting the steady-state probability Eq. and 
the transition rate matrix 1(61(1 to Eq. 1)560. one finds 



-r r 



G 2 (z,t) 



1 



4 h{z) 



(l + f 2 (z)) 2 e +r ° T ^ 



-(1- / 2 W) 2 e- r ° r/2(z) 



(62) 



w here f 2 (z) = (Ts/T )Vz 2 + A 2 , with T s = ^T 2 - <5 2 = 
VT+T^and A 2 = S 2 /T 2 . 

Now, it is straightforward to calculate the cumulants. 
In the long time limit r 3> (7 for instance, in terms 
of the average current I 2 = 67+7^/(7+ + 7") and the 
Fano factor F 2 in 1)35(1 . the three lowest cumulants are 
given by 



«i(r) = I 2 t, k 2 (t) = eI 2 F 2 r, 
ks(t) = e 2 I 2 (3{F 2 - 1/2) 2 + l/4)r, 



(63) 



where the electron charge (— e) is revived. These are in 
agreement with the phase-coherent quantum-mechanical 
results [lQ]. It is convenient to discuss the asymmetry 
(skewness) by the ratio A = n 3 /e 2 Ki (r — > 00), noticing 
the Fano factor F = lim^a, k 2 (t) / e«i (r) . The factor 
A 2 = 3(F 2 — 1/2) 2 + 1/4 is positive definite (positive 
skewness) and bounded by 1/4 < A 2 < 1. It is interesting 
to notice that A 2 is a monotonic function of F 2 and has 
its minimum A 2 — 1/4 for the minimum F 2 = 1/2 and 
maximum A 2 = 1 for the maximum F 2 = 1. Notice A = 
1 for a Poissonian and A = for a Gaussian distribution. 
Therefore, the dependence of A 2 on the gate charge N G 
is similar to that of the Fano factor F 2 (see Fig. with 
dips at N$ p in Eq. (gSJ. Together with Eq. it 
implies that the effective shot noise and the asymmetry 
of the probability distribution per unit charge transfer 
have their respective minimum values at the gate charge 
Nq P which depends on the tunnel-junction asymmetry 
and the interaction strength of the leads. 




FIG. 9: Contour of Eq. j57fr . The arguments of f 2 (z) = 
(Fs/F Q )Vz 2 + A 2 are tt/2, -tt/2, tt/2, -7r/2 along the branch 
cuts Ci,C2,C3, and Ca, respectively. 



The integral in Eq. I(57|) is along the contour depicted 
in Fig. |21 Notice that the contributions from the part 
along the branch cuts are zero and we are left with the 
multiple poles at z = 0. By residue theorem, the two- 
state probability P 2 (M, r) is given by Eq. (|BTjj) . 

The exact expression of P 2 (M, r) is cumbersome. For 
symmetric tunneling barriers with Ng = 1/2 (S = 
0, 7 + = 7~ = r ), however, Eq. ((fj2"|l reduces to 



G ( 2 s) (z,r) 



-Tot 

Az 



((l + z) 2 e 



-{l-zfe 



-t tz 



(64) 



Accordingly, P^' (M, r) is concisely given by 

"1 (ror) 2 ^- 1 , (r r) 2M 



F 2 (s) (M,r) = e- TaT 



2 (2M- 1)! 

1 (r r) 2Af + 1 

2 (2M+ 1)! 



(2M)! 
, for M > 1, (65) 



with P 2 (s) (0,r) = e- r ° T (l + r r/2), in agreement with 
Eq. (24) of Ref. [73 While this distribution resembles 
a sum of three Poisson distributions, it is not exactly 
Poissonian. 

For a highly asymmetric junctions R 3> 1 (7+ ^> 7 ~), 
the first term in Eq. I(62|l dominates the dynamics of 
G 2 (z,t) and its derivatives, and the characteristic func- 
tion is approximated by 

G ( 2 a) (z, t) « exp (-tT + rv/(7+7-)z 2 + 5 2 ) . (66) 

Now, the solution of P^ (M, r) is calculated by this 
equation and Eq. 1(60(1 . The leading order approxima- 
tion in 7~/7 + leads to the Poisson distribution, 



P 2 (a) (M,r) 



(7- 



Ml 



(67) 
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For a single tunneling-barrier, the charges are trans- 
ported by the Poisson process (F = 1) 45] . Therefore, we 
recover the Poisson distribution in the limit of strongly 
asymmetric junctions and in the regime of the two-state 
process, in which electrons see effectively single tunnel- 
barrier. 

For the intermediate barrier asymmetry, therefore, the 
probability Pi (M, r) of a two-state process is given by a 
distribution between Eq. (|65|l (for symmetric-junctions) 
and the Poissonian 1)6 7|) (for the most asymmetric- 
junctions). 



B. Four-state process; 2e p < eV < 2Ec 

Since we focus the FCS analysis on the case Nq = 1/2, 
the next simplest case to study is a four-state-model 
rather than the three-state-model discussed in the con- 
nection of the shot noise in the previous section. 

In the bias regime where the transport is governed by 
Coulomb blockade (eV < 2Ec) and yet the plasmons 
play an important role (eV > 2e p ), it is a fairly good ap- 
proximation to include only the four states with N — 0,1 
and n\ = 0, 1 (n rn = for m > 2). For general asymmet- 
ric cases, the rate matrix T(9) in Eq. H54|) is determined 
by ten participating transition rates (four rates from each 
junction, and two relaxation rates). The resulting eigen- 
values of \g(0, t)) are the solutions of a quartic equation, 
which is in general very laborious to solve analytically. 

For symmetric junctions (R — 1) with no plasmon re- 
laxation, however, the rate matrix is simplified to 



Too + 7io 

701 + 7n 

-7oo -701 

-710 -711 



with the matrix elements 



T R (0,t «- l,j). 
given by 



p 4 (s) M)h 5 (z,t = o)) = 



7ij 



-7oo^ 

7oo + 7io 


= r L (i, 



-701Z 

-7n^ 2 


7oi + 7n 



0,3) 



(68) 



The steady-state probability is then 



2(7oi +710) 



7oi 
7io 
701 
7io 



(69) 



Solving Eq. (|56|l with this probability and the rate ma- 
trix l(6"%)) is laborious but straightforward and one finds 



Q^ z T ^ _ g-(7oo+7lo+7oi+7ll W2 

,(! + *) 



{GlM ii±_EL } + { __ 4 



(70) 



with 



GAZ,T) = ( e (700,+7ll^+/4( Z ))r/2 ( l + Mz)\ 

I V U( Z )J 

A(z) - h(z) 



x 1- 



2(701 + 710)2 



+ <h(z)^-U(z)l (71) 



where A(z) and fi{z) are given by 

A(z) = 700 + 710 - 7oi - 7n - (7oo - 7n - 2701)2: 
h{z) = V(7oo - 7n) 2 - 4710701 V(z - a) 2 + b 2 (72) 

with dimensionless parameters a, b given by 

_ (700 - 7n) (7oo + 7io - 7oi - 7n) 



b = 



(700 - 7n) 2 - 47io7oi 
V47io7oi(7oq + 710 - 701 - 711)) 
(700 - 7n ) 2 - 4 7io7oi 



Integral of Gj{z, r) along the contour \z\ — 1 contains two 
branch points at z c — a ± ib, however, the integral along 
the branch cuts cancel out due to the symmetry under 
[fi(z) — > — fi(z)\. Therefore, the contribution from the 
branch cuts due to Gi(z, r) and Gi(—z, r) is zero to the 
probability (M, r) , and it is given by the residues 
only at z = 0, i.e. by Eq. I|60|l . The explicit expression 

of Pj s ' > (M, r) is cumbersome. 

The probability distribution P2 (M, r) for the two-state 
process deviates from Eq. I|65|l as a function of the asym- 
metry parameter R and reaches Poissonian in the case 
of strongly asymmetric-junctions. In a similar manner, 
P 4 (s) (M, t) deviates from P 2 (s) (M, r) as a function of ratio 



of the transition rates 



Hi- 



C. Numerical results 

It is worth mentioning that for strongly asymmetric 
junctions P(M,t) is Poissonian in the very low bias 
regime (eV < 2e p ), as seen in Eq. l|(77|) . It exhibits a 
crossover at eV = 2e p : P(M, r) deviates from Poisson 
distribution for 2e p < eV < 2Ec while it is Poissonian 
for eV < 2e p (at T = 0), as shown by the shot noise 
calculation. 



1. Voltage dependence 

The analytic results presented above are useful in in- 
terpreting the numerical results in Fig. I1UI where prob- 
ability P(M, t) for (a) symmetric junctions (R = 1) and 
(b) R = 100, in the case of LL parameter g — 0.5 with 
no plasmon relaxation (7^ = 0), is shown as a function 
of eV and M that is the number of transported electrons 
to the right lead during r such that during this time 
(M) = I c t = 10 electrons have passed to the right lead 
at eV = 2E C . 

The peak position of the distribution of P(M, r) is 
roughly linearly proportional to the average particle flow 
(M), and the width is proportional to the shot noise but 
in a nonlinear manner. In a rough estimate, therefore, the 
ratio of the peak width to the peak position is propor- 
tional to the Fano factor. Two features are shown in the 
figure. First, the average particle flow (the peak position) 
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the particle current at eV = 2Ec with no plasmon relaxation 
(7 P = 0); (a) for symmetric junctions (R = 1) and (b) for 
a highly asymmetric junctions (R — 100). Here g — 0.5, 
Ng = 1/2, and T = 0. Inset shows cross-sectional image of 
P(M, t) (blue solid line) as a function of M and the reference 
distribution function (a) Eq. 16511 (magenta dashed line) and 
(b) the Poisson distribution Eq. 1671 (red dashed), at (I) 
eV = 0.5E C , (II) eV = 1.5E a , and (III) eV = 2E C . 



runs with different slope when the bias voltage crosses 
new energy levels, i.e. at eV = 2e p and eV — 2Ec, that 
is consistent with the I — V study (compare Fig. llUf b) 
with Fig. EH). Notice E c = e p /g = 2e p for g = 0.5. Sec- 
ond, the width of the distribution increases with increas- 
ing voltage, with different characteristics categorized by 
eV = 2e p and eV = 2Ec- Especially in the bias regime 
eV > 2Ec in which the charge fluctuations participate 
to the charge transport, for the highly asymmetric case, 
the peak runs very fast while its width does not show 
noticeable increase. It causes the dramatic peak struc- 
ture in the Fano factor around eV — 2Eq as discussed in 



section Ivl 

The deviation of the distribution of probability 
P(M, r) due to the non-equilibrium plasmons from its 
low voltage (equilibrium) counterpart is shown in the in- 
sets. Notice in the low bias regime eV < 2e p , it follows 
Eq. (|65|l for the symmetric case (llOf a). inset I), and the 
Poissonian distribution (|67l) for the highly asymmetric 
case (|l(~)f b). inset I). The deviation is already noticeable 
at eV = 3e p = l.5E c for R = 1 (inset II inlOTa)). while 
it deviates strongly around eV = 2Ec for R = 100 (inset 
III inlMa)). 



2. Interaction strength dependence 



P(M,T) 



0.2 n 




FIG. 11: (color online) The probability P(M,t) that M elec- 
trons have passed through the right junction during the time 
r = 10//o, where Jo is the particle current with no plas- 
mon relaxation (7 P = 0); (a) with no plasmon relaxation 
(7 P = 0) and (b) with fast plasmon relaxation (7 P = 10 4 ). 
Here eV = 2E C , R = 100, iV G = 1/2 and T = 0. 

We have concluded in section that shot noise shows 
most dramatic behavior around eV = 2Eq due to in- 
terplay between the non-equilibrium plasmons and the 
charge fluctuations. To see its consequence in FCS, we 
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plot in Fig. 1111 the probability P(M, r) as a function of 
the particle number M and the interaction parameter g 
for t = 10e/I(eV = 2Ec,j P = 0) (a) with no plasmon 
relaxation and (b) with fast plasmon relaxation. 

The main message of Fig. Illf a) is that the shot noise 
enhancement, i.e. the broadening of the distribution 
curve, is significant in the strong interaction regime with 
gradual increase with decreasing g. Fast plasmon relax- 
ation consequently suppresses the average current dra- 
matically as shown in Fig. II If implying the Fano factor 
enhancement is lost. Effectively, the probability distribu- 
tion of P(M, t) for different interaction parameters maps 
on each other almost identically if the time duration is 
chosen such that (M) equals for all g. 

VII. CONCLUSIONS 

We have studied different transport properties of a 
Luttinger-liquid single-electron transistor including av- 
erage current, shot noise, and full counting statistics, 
within the conventional sequential tunneling approach. 

At finite bias voltages, the occupation probabilities of 
the many-body states on the central segment is found 
to follow a highly non-equilibrium distribution. The en- 
ergy is transferred between the leads and the quantum 
dot by the tunneling electrons, and the electronic iden- 
tity is dispersed into the plasmonic collective excitations 
after the tunneling event. In the case of nearly symmet- 
ric barriers, the distribution of the occupation probabil- 
ities of the non-equilibrium plasmons shows impressive 
contrast depending on the interaction strength: In the 
weakly interacting regime, it is a complicated function 
of the many-body occupation configuration, while in the 
strongly interacting regime, the occupation probabilities 
are determined almost entirely by the state energies and 
the bias voltage, and follow a universal distribution re- 
sembling Gibbs (equilibrium) distribution. This feature 
in the strong interaction regime fades out with the in- 
creasing asymmetry of the tunnel-barriers. 

We have studied the consequences of these non- 
equilibrium plasmons on the average current, shot noise, 
and counting statistics. Most importantly, we find that 
the average current is increased, shot noise is enhanced 
beyond the Poisson limit, and full counting statistics 
deviates strongly from the Poisson distribution. These 
non-equilibrium effects are pronounced especially in the 
strong interaction regime, i.e. g < 0.5. The overall trans- 



port properties are determined by a balance between phe- 
nomena associated with non-equilibrium plasmon distri- 
bution that tend to increase noise, and charge fluctu- 
ations that tend to decrease noise. The result of this 
competition is, for instance, a non-monotonic voltage de- 
pendence of the Fano factor. 

At the lowest voltages when charge can be transported 
through the system, the plasmon excitations are sup- 
pressed, and the Fano factor is determined by charge 
fluctuations. Charge fluctuations are maximized when 
the tunncling-in and tunncling-out rates are equal, which 
for symmetric junctions occurs at gate charge Nq = 1/2. 
At these gate charges the Fano factor acquires its lowest 
value which at low voltages is given by a half of the Pois- 
son value, known as 1/2 suppression, as only two states 
are involved in the transport, at somewhat larger voltages 
increases beyond the Poisson limit as plasmon excitations 
are allowed, and at even higher voltages exhibits a local 
minimum when additional charge states are important. 
If the non-equilibrium plasmon effects are suppressed e.g. 
by fast plasmon relaxation, only the charge fluctuation 
effects survive, and the Fano factor is reduced below its 
low-voltage value. The non-equilibrium plasmon effects 
are also suppressed in the non- interacting limit. 
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APPENDIX A: THE TRANSITION 
AMPLITUDES IN THE QUANTUM DOT 

In this appendix, we derive the transition amplitudes, 
Eqs. i|15|) and (JTSJ. As shown in Eq. Q14[l . the zero-mode 
overlap of the QD transition amplitude is unity or zero. 
Therefore, we focus on the overlap of the plasmon states. 
It is enough to consider 

|({n'}|vI/t D (^)|{n})| 2 , X t = X L ,X R 

due to the symmetry between matrix elements of 
tunneling-in and tunneling-out transitions, as in Eq. 

dm, 



|({n'}|vI/t(A^)|{n})| 2 = [l(K}|^(^)IW)| 2 + <{n'}|^(X,)|{n})(W|^_ r (X,)|K}) 



= 2\{{n'M{X e )\{n})f 



(Al) 



where r = +(— ) denotes the right(left)-moving compo- nents cancel out due to fermionic anti-commutation re- 
nent, and the cross terms of oppositely moving compo- 
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lations. 

The transition amplitudes at is identical to that 
at Xr. For simplicity, we consider the case at Xl = 
only. The overlap elements of the many-body oc- 
cupations ({n}\ — (rii, ri2, • • • , n m , ■ ■ ■ | and = 
\ n 'i, n 2, ■■■ , n 'm>'--) are 

1 OO 

\({n}\4(x = 0)|{n'})| 2 = — J] \(n m \^ m \n'J\ 2 , 



(A2) 
is 



where ip m = exp [\ m (b m + 6}„)] with X rn = —-——^ 
the bosonized field operator at an edge of the wire with 



open boundary conditions (see for instance Ref. |37 
Here g is the interaction parameter, m is the mode index 
(and the integer momentum of it) , and M is the number 
of transport sectors; if M > 1, the contributions of the 
different sectors must be multiplied. The operators b m 
and b\ n denote plasmon annihilation and creation and A 
is a high energy cut-off. 

Using the Baker-Haussdorf formula 



ip m = exp [A m (a m + a*,)] 
and the harmonic oscillator states 

|„) = LJ f |0>, (n| = 
v nl 

one can show that, if n < n', 

U|2(n-n') 



(A3) 



^ (A4) 



nl 



| (n\cp\n') | 2 = e |A| 



n\n'l \(n — n')\ y 
x $(?i + l,n-n' + l;-|A| 2 ) 5 



(A5) 



where $(2:, a/; z) is a degenerate hypergeometric function 
defined by [Hi 



$(a;, x';z) 



U=0 



z £ (a;-l+£)! (x'-l)! 
J\ (ar-l)! (a;' -!+£)! 



(A6) 



If n' < n, the indices n and n' are exchanged in Eq. 1A5I . 
The function (f>(x,a;'; z) is a solution of the equation 

z<9 2 $ + (x - z)d z <f> - x$ = 0. 

By solving this differential equation with the proper nor- 
malization constant, one obtains 

<D(n+l,n-n'+l;-|A| 2 ) - e~W' L^ n ' (\X\ 2 ), 

where (y) is the Laguerre polynomials |7l| . In terms of 
the Laguerre polynomials, therefore, the transition am- 
plitude (| A5|) is written by 



\(n m \<p m \n' m )\ 2 = 



-1/gmM 



,!<)! 



(pM)K-»J n (»! 



\n' m 



\ gmM 



n 2 



(A7) 



where = min(n',n) and n^>) — max(n',n). 

We introduce a high frequency cut-off m c ~ kpLn/Tt 
to cure the vanishing contribution due to e _1 / 9mM , 



Trie 



/4mg 



2?rA -"-I " 2i D 

m— 1 



where the exponent is a = (g 1 — 1)/M. We arrive at the 
desired form of the on-dot transition matrix elements, 



\{n'}\^ a (X e )\{n})\ 2 
1 



1 

To 



7rA 
1^ 



nn 



g u mM 



W m -n m \ «), 



L 



1 



g v mM 



(A9) 



APPENDIX B: UNIVERSAL OCCUPATION 
PROBABILITY 

In this appendix, we derive the universal distribution 
of the occupation probability Eq. (j26|) . 

Since the occupation probability of the plasmon many- 
body states is a function of the state energy in the 
strong interaction regime, we introduce the dimension- 
less energy n = ^ m ■ n rn of the state with {n} — 
(m, ri2, • • • , n m , • • • ) plasmon occupations. Excluding 
the zero mode energy, therefore, the energy of the state 
{n} is given by Ejj({n}) — Ejj(n) = ne p with the state 
degeneracy D(n), i.e. the number of many-body states 
{n} satisfying n — m ■ n m , asymptotically following 
the Hardy-Ramanujan formula 



D(n) 



/2n/3 



asyr 

m 

/(4VSr, 



(Bl) 



We denote n sd by the corresponding dimcnsionlcss bias 
voltage eV — n sd e p . 

We obtain an analytic approximation to the occupa- 
tion probability P(n) at zero temperature by setting the 
on-dot transition elements in l|15f) to unity and consid- 
ering the scattering-in and scattering-out processes for a 
particular many-body state {n}. 

The total scattering rates at zero temperature are 
given by a simple power-law Eq. (1131) , 

r(n <— m) — 0(— rt + m — n sd /2)(—n + m — n sd /2) a 
+ 6(-n + m + n sd /2)(-n + m + n sd /2) a 
« Q(-n + m - n sd /2)(-n + m- n sd /2) a (B2) 

where the constant factor in Eq. I|13|) is set to unity. 
The master equation now reads 

d 

— P(n) =^2[P(m)D(m)T(m -> n)-P(n)D(m)T(n -> m 



(B3) 



20 



To solve this master equation, we assume an ansatz of a 
power-law 



P(m) = P(n)q« 



(B4) 



In the steady-state, master equation (|B3|) in terms of this 
ansatz becomes 



(m-n + n sd /2) a D(m)q™ 



n+n Bd /2 

« J2 (n-m' + n sd /2) a D{m'), (B5) 

m'=0 

in which the sum in the LHS runs from rrii = max(0, n — 
n sd /2), where maxfi, x') returns larger value of x and x' . 
Using a saddle point integral approximation 

if/'(xo) = 0, f'(x o )<0, (B6) 

we solve Eq. (|B5|) to obtain an equation for lng„ and 
find that for a large n, 



where 



exp(z)/z w \/n exp(F(n)), 

"7 L^a — 



(B7) 



and P(n) is a slowly varying function of n for n>l, 



F(n 



2 7T 



hi 



V6 



1 



3 C a ^Jn 

We assume an ansatz for the solution of z in Eq. (|B7J| . 



n \4 a + 1 

(B8) 



= (hm)/2 + F(n) + ^/(lnn)^ + F(n) - K + i] (B9) 



and find a constant F which minimizes the correction 
term r\. Putting this ansatz into Eq. (|B7|I and solving 
the equation for rj, we find at K w 0.8, the correction 
term 77 is negligibly small (77 rs 0.01). 



Noting /(n) = \/ (In n) /2 + F(n) — if is almost linear 
function in the regime of our interest (3 < n < 15), we 
linearize it around a value of interest n — hq (for instance, 
n = 9), 

/(n) w f'(n )(n - n Q ) + f(n ), 



and solve Eq. I|B7|) by ansatz (|B9|) with above linearized 
form; 



Cq. 2 



In n + F(n) + f'(n )n + /(no) - /'(n )n . 

(BIO) 

Apply <9 m ln[P(m)] « ln[g„] to Eq. (|B10|i . and solve 
the integral equation for ln[P(n)], 



ln[P(n)] « -C Q 



(lnn)/2 + F(n) 



+ /'(no)n + /(no)-no/'M • (Bll) 



The leading order approximation p(°) (n) to the probabil- 
ity P U niv (n) of the average occupation from this integral 
results in Eq. (2| 



P(°)(n) = Z- 1 rT 1 ^" = Z- 1 • e -i^ nlogn , (B12) 



where Z is the partition function. 

A more accurate approximation PW(n) can be derived 
by solving the integral Eq. (|B11|) to a higher degree of 
precision, which yields 



(»+=?)) 



-CW2 / ^6 

x — C Q 



-C a n 



exp —C a n 



F(n )-K- 



ln3 



x exp 



/ 4tt _\ (n-n sd /2\ 2n/n ° d / n„,\ / 2 /n sd \ 



(B13) 



where no = 9 is used and f{rio) w f(no) 2 with minor 
correction is utilized for formal simplicity. Note the first 
term with normalization constant Z approaches P^ (n) 
in Eq. (|B12() as n/n sd — ► 0, noticing C Q = 2 (a + l)/n s d- 



The integral equation (|B 1 1|> can be solved even without 
approximating on f(n), with the expense of more cum- 
bersome appearance of P(n). 
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